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Amendments to the Claims: 

Please amend claims 34 and 71 as provided below. 

This listing of claims will replace all prior versions, and listings, of claims in the 
application: 

Listing of Claims: 

1 . (Previously Presented) A computer readable medium or media having stored thereon 
computer-implemented instructions causing a processor to perform the steps of: 

(a) providing a data structure relating a plurality of reactants to a plurality of 
reactions of a biochemical reaction network, wherein each of said reactions 
comprises a reactant identified as a substrate of the reaction, a reactant identified 
as a product of the reaction and a stoichiometric coefficient relating said substrate 
and said product, and wherein at least one of said reactions is a regulated reaction; 

(b) providing a constraint set for said plurality of reactions, wherein said constraint 
set comprises a variable constraint for said regulated reaction; 

(c) determining at least one flux distribution that minimizes or maximizes an 
objective function when said constraint set is applied to said data structure, 
wherein said at least one flux distribution determines a systemic property of said 
biochemical reaction network, and wherein said systemic property is dependent 
upon the flux through said regulated reaction; and 

(d) providing information resulting from steps (a) through (c) to a user. 

2. (Original) The computer readable medium or media of claim 1 , wherein said variable 
constraint is dependent upon the outcome of at least one reaction in said data structure. 

3. (Original) The computer readable medium or media of claim 1 , wherein said variable 
constraint is dependent upon the outcome of a regulatory event. 

4. (Original) The computer readable medium or media of claim 1, wherein said variable 
constraint is dependent upon time. 
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5. (Original) The computer readable medium or media of claim 1, wherein said variable 
constraint is dependent upon the presence of a biochemical reaction network participant. 

6. (Original) The computer readable medium or media of claim 5, wherein said participant 
is selected from the group consisting of a substrate, product, reaction, protein, 
macromolecule, enzyme and gene. 

7. (Original) The computer readable medium or media of claim 1, wherein said 
biochemical reaction network comprises metabolic reactions. 

8. (Original) The computer readable medium or media of claim 1 , further comprising a 
regulatory data structure, wherein said variable constraint is dependent upon an outcome 
of a regulatory event represented by said regulatory data structure. 

9. (Original) The computer readable medium or media of claim 8, wherein said regulatory 
data structure represents a regulatory event selected from the group consisting of 
transcription of a gene, translation of an RNA, post-translational modification of a 
protein, inhibition of a protein, activation of a protein, assembly of a protein, change in 
pH, change in redox potential, change in temperature, passage of time, and degradation of 
a protein. 

1 0. (Original) The computer readable medium or media of claim 8, wherein said regulatory 
event is due to a signal transduction pathway. 

1 1 . (Original) The computer readable medium or media of claim 8, wherein said 
biochemical reaction network and said regulatory data structure represent reactions or 
events that occur in a single cell. 

12. (Original) The computer readable medium or media of claim 8, wherein said 
biochemical reaction network represents reactions that occur in a first cell in a population 
of cells and said regulatory data structure represents events that occur in a second cell in 
said population. 
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13. (Original) The computer readable medium or media of claim 12, wherein said population 
of cells comprises cells of a multicellular organism. 

14. (Original) The computer readable medium or media of claim 1, further comprising a 
constraint function that correlates an outcome of a regulatory event with said variable 
constraint. 

15. (Original) The computer readable medium or media of claim 14, wherein said constraint 
function is binary. 

16. (Original) The computer readable medium or media of claim 14, wherein said regulatory 
event is represented by Boolean logic. 

17. (Canceled) 

18. (Previously Presented) The computer readable medium or media of claim 1 , further 
comprising a range of said feasible flux distributions. 

19. (Previously Presented) The computer readable medium or media of claim 1 , wherein said 
commands comprise an optimization problem. 

20. (Original) The computer readable medium or media of claim 1 9, wherein said 
optimization problem comprises a linear optimization problem or a nonlinear 
optimization problem. 

21 . (Previously Presented) The computer readable medium or media of claim 1 , further 
comprising a user interface capable of sending at least one command for modifying said 
data structure, said constraint set or said commands for applying said constraint set to 
said data representation, or a combination thereof. 

22. (Original) The computer readable medium or media of claim 2 1 , wherein said user 
interface further comprises links which a user may select to access additional information 
relating to said plurality of reactions. 
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23. 
24. 
25. 
26. 

27. 
28. 

29. 
30. 
31. 

32. 



(Original) The computer readable medium or media of claim 1, wherein said data 
structure comprises a set of linear algebraic equations. 

(Original) The computer readable medium or media of claim 1, wherein said data 
structure comprises a matrix. 

(Original) The computer readable medium or media of claim 1, further comprising 
commands for representing said at least one flux distribution as a flux distribution map. 

(Original) The computer readable medium or media of claim 1, wherein at least one 
reactant in said plurality of reactants or at least one reaction in said plurality of reactions 
is annotated. 

(Original) The computer readable medium or media of claim 26, wherein said annotation 
comprises assignment of said at least one reactant to a compartment. 

(Original) The computer readable medium or media of claim 27, wherein a first 
substrate or product in said plurality of reactions is assigned to a first compartment and a 
second substrate or product in said plurality of reactions is assigned to a second 
compartment. 

(Original) The computer readable medium or media of claim 26, wherein said annotation 
comprises assignment to an open reading frame or protein. 

(Previously Presented) The computer readable medium or media of claim 26, wherein 
said annotation comprises a confidence limit for occurrence of a reaction. 

(Original) The computer readable medium or media of claim 1, further comprising a 
gene database relating one or more reactions in said data structure with one or more 
genes or proteins in particular organism. 

(Original) The computer readable medium or media of claim 1, wherein said 
biochemical reaction network comprises reactions that are selected from the group 
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consisting of glycolysis, the TCA cycle, the pentose phosphate pathway, respiration, 
biosynthesis of an amino acid, degradation of an amino acid, biosynthesis of a purine, 
biosynthesis of a pyrimidine, biosynthesis of a lipid, metabolism of a fatty acid, 
biosynthesis of a cofactor, metabolism of a cell wall component, transport of a metabolite 
and metabolism of carbon, nitrogen, sulfur, phosphate, hydrogen or oxygen. 

33. (Original) The computer readable medium or media of claim 1 , wherein a plurality of 
said reactions are regulated reactions and said constraints for said regulated reactions 
comprise variable constraints. 

34. (Currently Amended) A method for determining a systemic property of a biochemical 
reaction network, comprising: 

(a) providing a data structure relating a plurality of reactants to a plurality of 
reactions of a biochemical reaction network, 

wherein each of said reactions comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the reaction and a stoichiometric coefficient 
relating said substrate and said product, and wherein at least one of said reactions is a 
regulated reaction; 

(b) providing a constraint set for said plurality of reactions, 

wherein said constraint set comprises a variable constraint for said regulated reaction; 

(c) providing a condition-dependent value to said variable constraint; 

(d) providing an objective function; 

(e) determining on a computer at least one flux distribution that minimizes or 
maximizes said objective function when said constraint set is applied to said data 
structure, 

wherein said at least one flux distribution is determinative of a systemic property of said 
biochemical reaction network; and 

(f) providing said systemic property of said biochemical reaction network to a user. 
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35. 

36. 

37. 
38. 

39. 
40. 
41. 

42. 



(Previously Presented) The method of claim 34, wherein said condition-dependent value 
provided to said variable constraint comprises a value conditioned on the outcome of at 
least one reaction in said data structure. 

(Previously Presented) The method of claim 34, wherein said condition-dependent value 
provided to said variable constraint comprises a value conditioned on the outcome of a 
regulatory event. 

(Previously Presented) The method of claim 34, wherein said condition-dependent value 
provided to said variable constraint comprises a value conditioned on time. 

(Previously Presented) The method of claim 34, wherein said condition-dependent value 
provided to said variable constraint comprises a value conditioned on the presence of a 
biochemical reaction network participant. 

(Original) The method of claim 38, wherein said participant is selected from the group 
consisting of a substrate, product, reaction, en2yme, protein, macromolecule and gene. 

(Original) The method of claim 34, wherein said biochemical reaction network 
comprises metabolic reactions. 

(Previously Presented) The method of claim 34, wherein said data structure further 
comprises a regulatory data structure, wherein said value provided to said variable 
constraint is changed due to an outcome of a regulatory event represented by said 
regulatory data structure. 

(Original) The method of claim 41, wherein said regulatory event is selected from the 
group consisting of transcription of a gene, translation of an RNA, post-translational 
modification of a protein, inhibition of a protein, activation of a protein, assembly of a 
protein, change in pH, change in redox potential, change in temperature, passage of time, 
and degradation of a protein. 
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43. 
44. 
45. 
46. 

47. 

48. 

49. 
50. 

51. 

52. 
53. 



(Original) The method of claim 41, wherein said regulatory event is due to a signal 
transduction pathway. 

(Original) The method of claim 41, wherein said biochemical reaction network and said 
regulatory data structure represent reactions or events that occur in a single cell. 

(Original) The method of claim 41, wherein said regulatory event comprises a regulatory 
reaction. 

(Original) The method of claim 41, wherein said biochemical reaction network 
represents reactions that occur in a first cell in a population of cells and said regulatory 
data structure represents events that occur in a second cell in said population. 

(Original) The method of claim 46, wherein said population of cells comprises cells of a 
multicellular organism. 

(Original) The method of claim 41, further comprising a constraint function that 
correlates an outcome of a regulatory event with said variable constraint. 

(Original) The method of claim 48, wherein said constraint function is binary. 

(Original) The method of claim 48, wherein said regulatory event is represented by 
Boolean logic. 

(Original) The method of claim 48, wherein said constraint function correlates a first set 
of outcomes of said regulatory data structure with a first binary value and a second set of 
outcomes of said regulatory data structure with a second binary value. 

(Original) The method of claim 48, wherein said constraint function correlates a set of 
outcomes of said regulatory data structure with a single integer value. 

(Original) The method of claim 34, wherein said flux distribution is determined by 
optimization. 
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54. 

55. 

56. 

57. 
58. 

59. 
60. 



(Original) The method of claim 53, wherein said optimization comprises linear 
optimization or non linear optimization. 

(Original) The method of claim 34, further comprising a step of modifying said data 
structure or said constraint set, or both. 

(Original) The method of claim 34, wherein said data structure comprises a set of linear 
algebraic equations. 

(Original) The method of claim 34, wherein said data structure comprises a matrix. 

(Original) The method of claim 34, further comprising a step of producing a flux 
distribution map. 

(Original) The method of claim 34, wherein said biochemical reaction network 
comprises reactions that are selected from the group consisting of glycolysis, the TCA 
cycle, pentose phosphate pathway, respiration, biosynthesis of an amino acid, degradation 
of an amino acid, biosynthesis of a purine, biosynthesis of a pyrimidine, biosynthesis of a 
lipid, metabolism of a fatty acid, biosynthesis of a cofactor, metabolism of a cell wall 
component, transport of a metabolite and metabolism of a carbon source, nitrogen source, 
oxygen source, phosphate source, hydrogen source or sulfur source. 

(Original) The method of claim 34, wherein said systemic property is selected from the 
group consisting of growth, energy production, redox equivalent production, biomass 
production, production of biomass precursors, production of a protein, production of an 
amino acid, production of a purine, production of a pyrimidine, production of a lipid, 
production of a fatty acid, production of a cofactor, production of a cell wall component, 
transport of a metabolite, development, intercellular signaling, and consumption of 
carbon nitrogen, sulfur, phosphate, hydrogen or oxygen. 

(Original) The method of claim 34, wherein said systemic property is selected from the 
group consisting of degradation of a protein, degradation of an amino acid, degradation 
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of a purine, degradation of a pyrimidine, degradation of a lipid, degradation of a fatty 
acid, degradation of a cofactor and degradation of a cell wall component. 

62. (Original) The method of claim 34, wherein said variable constraint comprises a 
condition-dependent constraint value and a constraint function, wherein said variable 
constraint is modified by said constraint function acting upon said condition-dependent 
constraint value. 

63. (Original) The method of claim 62, wherein said constraint function is binary. 

64. (Original) The method of claim 34, further comprising providing a gene database 
relating one or more reactions in said data structure with one or more open reading 
frames or proteins in a particular organism. 

65. (Original) The method of claim 64, further comprising identifying an open reading frame 
that encodes a protein that performs a reaction in said plurality of reactions. 

66. (Withdrawn) The method of claim 64, further comprising identifying a protein that 
performs a reaction in said plurality of reactions. 

67. (Withdrawn) A method for determining a phenotype of a mutant of an organism, 
comprising: 

(i) identifying a reaction that is not naturally present in a particular organism, and 

(ii) determining a systemic property of a biochemical reaction network according to 
the method of claim 34, wherein said data structure relates a plurality of reactants 
for said organism to a plurality of reactions of a biochemical reaction network of 
said organism and further comprises said reaction that is not naturally present in 
said organism. 

68. (Withdrawn) A method for determining a phenotype of a mutant of an organism, 
comprising: 
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(i) identifying a reaction that is related to an open reading frame or protein in said 
gene database, and 

(ii) determining a systemic property of a biochemical reaction network according to 
the method of claim 34, wherein said reaction that is related to said open reading 
frame or protein is not present in said data structure or is constrained to have no 
flux. 

69. (Withdrawn) A method for determining the effect of a drug on the activity of one or 
more reactions in a biochemical reaction network, comprising: 

(i) identifying a reaction that is related to an open reading frame or protein in said 
gene database; 

(ii) identifying a candidate drug that alters expression of said open reading frame or 
activity of said protein, and 

(iii) determining a systemic property of a biochemical reaction network according to 
the method of claim 34, wherein said reaction that is related to said open reading 
frame or protein is not present in said data structure, is constrained to have a 
reduced flux, or is constrained to have no flux. 

70. (Original) The method of claim 34, wherein a plurality of said reactions are regulated 
reactions and said constraints for said regulated reactions comprise variable boundary 
values. 

71 . (Currently Amended) A method for determining a systemic property of a biochemical 
reaction network at a first and second time, comprising: 

(a) providing a data structure relating a plurality of reactants to a plurality of 

reactions of a biochemical reaction network, 
wherein each of said reactions comprises a reactant identified as a substrate of the 
reaction, a reactant identified as a product of the reaction and a stoichiometric coefficient 
relating said substrate and said product, and wherein at least one of said reactions is a 
regulated reaction; 
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(b) providing a constraint set for said plurality of reactions, wherein said constraint 
set comprises a variable constraint for said regulated reaction; 

(c) providing a condition-dependent value to said variable constraint; 

(d) providing an objective function; 

(e) determining on a computer at least one flux distribution at a first time that 
minimizes or maximizes said objective function when said constraint set is 
applied to said data structure, thereby determining a systemic property of said 
biochemical reaction network at said first time; 

(f) modifying said value provided to said variable constraint; 

(g) repeating step (e) wherein said at least one flux distribution is determined at a 
second time, thereby determining a systemic property of said biochemical 
reaction network at a second time, and 

(h) providing said systemic property of said biochemical reaction network to a user at 
said first, second or first and second time. 

72. (Original) The method of claim 71 , wherein said value is modified based on said flux 
distribution at said first time. 

73. (Original) The method of claim 71 , wherein said value is modified based on a change in 
an environmental condition. 

74. (Original) The method of claim 71, further comprising repeating steps (e) through (g) for 
a specified number of timepoints. 
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